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Fluctuations of the atomic positions are at the core of a large class of unusual material 
properties ranging from quantum para-electricity to high temperature superconductivity. Their 
measurement in solids is the subject of an intense scientific debate focused on seeking a method¬ 
ology capable of establishing a direct link between the variance of the atomic displacements and 
experimentally measurable observables. Here we address this issue by means of non-equilibrium 
optical experiments performed in shot-noise limited regime. The variance of the time dependent 
atomic positions and momenta is directly mapped into the quantum fluctuations of the photon 
number of the scattered probing light. A fully quantum description of the non-linear interac¬ 
tion between photonic and phononic fields is benchmarked by unveiling the squeezing of thermal 
phonons in a-quartz. 

In a classical description the displacement of the atoms along the vibrational eigenmodes of a crystal can 
be measured with unlimited precision. Conversely, in the quantum formalism positions and momenta of the 
atoms can be determined simultaneously only within the boundary given by the Heisenberg uncertainty principle. 
For this reason, in real materials, in addition to the thermal disorder, the atomic displacements are subject to 
fluctuations which are intrinsic to their quantum nature. While various evidences suggest that such quantum 
fluctuations may be of relevance in determining the onset of intriguing material properties, such as quantum 
para-electricity, charge density waves, or even high temperature superconductivity [1-7], the possibility of 
measuring quantum fluctuations in solids is the subject of an intense debate [8-21]. 

The time evolution of atomic positions in materials is usually addressed by means of non-equilibrium optical 
spectroscopy. An ultrashort light pulse (the pump) impulsively perturbs the lattice and a second one (the probe), 
properly delayed in time, measures a response that is proportional to the spatially-averaged instantaneous atomic 
positions. In those experiments, the time dependent atomic displacements are often revealed by an oscillating 
response, commonly dubbed coherent phonon response [22-32], at frequencies characteristic of the vibrational 
modes of the material. In this framework, it has been shown that a non-linear light-matter interaction can 
prepare non classical vibrational states [8,13] such as squeezed states, where the fluctuations of the lattice 
position (or momentum) can fall below the thermal limit. A reduction below the vacuum limit is known as 
vacuum squeezing [16]. 

Here we propose a joint experimental and theoretical approach to access the fluctuations of the atomic 
positions in time domain studies. An experimental apparatus that allows for the measurement of the photon 
number quantum fluctuations of the scattered probe pulses in a pump and probe setup is adopted. The connection 
between the measured photon number uncertainty and the fluctuations of the atomic positions is given by a fully 
quantum mechanical theoretical description of the time domain process. Overall we prove that, under appropriate 
experimental conditions, the fluctuations of the lattice displacements can be directly linked to the photon number 
quantum fluctuations of the scattered probe pulses. Our methodology, that combines non-linear spectroscopic 
techniques with a quantum description of the electromagnetic fields, is benchmarked on the measurement of 
phonon squeezing in a-quartz. 
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Results 

Shot-noise limited pump and probe experiments 

In the non-linear spectroscopy formalism, the excitation mechanism of phonon states in transparent materials is 
called impulsive stimulated Raman scattering (ISRS) [29]. The susceptibility tensor connects the induced 
third order polarization to three fields; 82 ( 102 ) and 83 ( 103 ) [33]. In conventional two pulses pump 

and probe experiments, the Helds 82 ( 1 ^ 2 ) and £’ 3 (^ 3 ) are two different frequency components of the pump laser 
pulse. In particular, all photon pairs such that C 03 — L 02 = fi, where fl is the frequency of the Raman active 
vibrational mode, contribute to ISRS [34]. The interaction of the probe field 8 i(ioi) with the photo-excited 
material induces an emitted field, £’ef(<^) which depends on the pump-probe delay and carries information about 
the specific Raman mode excited in the crystal. Experimental details are reported in (Supplementary Note 
1). We choose a polarization layout designed to excite E-symmetry Raman optical modes in a-quartz at room 
temperature and get an emitted field with polarization orthogonal to the probe one [35] (Supplementary Fig. 1). 

We propose here a new approach to time domain studies. The experimental layout is similar to standard 
pump and probe experiments. The sample is excited by an ultrashort pump pulse and the time evolution of 
the response is measured by means of a second much weaker probe pulse, that interacts with the photo-excited 
material at a delay time r. The unique characteristics of our setup are; i) unlike standard experiments, where 
the response is integrated over many repeated measurements, our system can measure individual pulses; ii) the 
apparatus operates in low noise conditions allowing for the measurement of intrinsic photon number quantum 
fluctuations. In detail, we adopt a differential acquisition scheme where each probe pulse is referenced with a 
second pulse which has not interacted with the sample. For each measurement the differential voltage is digitized 
and integrated, giving the transmittance AT^ for the ***’ measurement. For every given pump and probe delay 
T, we repeat this single pulse measurements for N = 4000 consecutive pulses. Fig. 1 (a) gives a useful visual 
representation of the obtained data. For one pump and probe scan I the normalized histogram of N = 4000 
acquired pulses for each delay time is shown. Each histogram represents the distribution of the measured ATi for 
a specific delay time t. For a clearer visualization of the physically meaningful information in the time evolution 
of the statistical distribution. Fig. 1 (b) reports the histogram centered at zero. 

The pump and probe scan is repeated several times and each 8 ^ scan provides ATi, and 

AT^r = 'k [^Ti — Finally the averages of these two quantities are calculated over all M scans as 

AT — 4- AT^*) and AT — 4- AT^^'> 

mean — Z^/=l mean “UU var — Z^l=l var' 

The time domain response, averaged over M = 10 scans, is shown in Fig. 2 (a) for a representative pump 
fiuence of 14mJcm“^ (a pump fiuence dependent study is reported later). The blue curve depicts the time 
evolution of the mean value of the transmittance ATmean, whereas the red curve shows the time evolution of its 
variance ATvar. The Fourier transform of the mean (Fig. 2 (c), blue curve) has a single peak which is ascribed 
to the E-symmetry quartz vibrational mode with frequency fl = 128cm“^ = 3.84THz [36|. The same frequency 
component is observed in the Fourier transform of the variance (Fig. 2 (c), red curve). In addition, a second 
peak at twice the phonon frequency appears exclusively in the variance. A wavelet analysis of the variance 
oscillations allows for a time domain study of the two frequency components (Fig. 2 (b)); one notices that, 
while the fundamental frequency survives for roughly 7 ps, the 281 component vanishes within the first 2ps. The 
different lifetimes between the fl and 2 fl components of the variance are seen also by a close inspection of the 
raw data distribution plotted in Fig. 2 (b). 

Note that the 2fl in our data is visible only in experimental conditions where the noise is dominated by 
the quantum uncertainty, a situation which is known as shot-noise regime. In such conditions ATvar measures 
the quantum variance of the scattered probe photon number. A full characterization of the detection system 
is reported in (Supplementary Note 2), including the shot-noise characterization (Supplementary Fig. 2) and 
the analysis of classical noise sources (Supplementary Fig. 3 and Supplementary Fig. 4). It should further be 
stressed that in experimental conditions where the noise is larger and dominated by classical sources the 2 fl 
contribution to the noise becomes unmeasurable. 

The presence of the 2fl frequency component is suggestive of phonon squeezing, as it has been indicated by 
Raman tensor models [8,9,12[. Nevertheless, the experimental evidences up to date lack a direct comparison 
with a reliable quantum noise reference [10,11,13, 20[. Hence, in these experiments the observation of the 2fl 
frequency in the optical noise is considered as an indication of phonon squeezing, but not an unequivocal proof. 
In details, a 2fl oscillating optical noise was reported in [13[, but later ascribed to an artifact [15[ due to the 
experimental amplification process. In particular, it has been demonstrated that amplification artifacts become 
more relevant when, using a lock-in amplifier based acquisition, the time constant of the lock-in increases with 
respect to the time between steps in the pump-probe delay. This gives rise to maxima in the noise where the 
derivative of the mean signal is maximal [15[. Here we use a pump power density which is almost 3 orders of 
magnitude higher than in [13,15[. In addition, we observe a 2fl frequency component in the optical variance 
which exhibits maxima in correspondence with the minima of the derivative of the mean signal, hence ruling out 
possible artifacts of the kind described in [15[. 
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Fully quantum description of impulsive stimulated Raman scattering 

In order to predict how the fluctuations of the atomic positions in a lattice can be mapped onto the photon 
number quantum fluctuations of the probe field, we develop a novel theoretical approach to time domain studies 
which treats quantum mechanically both the material and the optical fields involved in the non-linear processes. 
Several semiclassical models describe the possibility of generating "classical" (coherent states) and non classical 
vibrational states by photo-excitation. In particular, for transparent materials like quartz, the most commonly 
used approach is to adopt Raman tensor models where the interaction between photons and phonons is not 
mediated by dipole allowed electronic transitions. In this condition, interactions linear in the phonon operators 
allow for the generation of coherent vibrational states, while high order interactions are required for the generation 
of non classical squeezed states [9,12,32]. In materials with allowed dipole transitions, as in presence of excitons, 
different models based on electron-phonon coupling Hamiltonians have been proposed. In those models it has 
been shown that squeezed phonon states can result only by successive excitations with a pair of pulses [16,17]. 
All these models mainly adopt semiclassical approaches where the optical fields are described classically [33], and 
therefore are unable to reproduce the quantum proprieties of the probe optical field that can be measured with 
the shot-noise limited pump and probe setup presented here. The key aspect of our approach, allowing us to 
bridge this gap, is to study both generation and detection of phonon states using a fully quantum formalism 
through an effective photon-phonon interaction, which is descriptive of experiments in transparent systems, such 
as a-quartz. The basic tool is a quantum Hamiltonian able to describe both pump and probe processes. Being 
linear and bilinear in the photon and phonon operators, this Hamiltonian accounts for the possible generation 
of coherent and squeezed phonon states through the pump process. In particular, it models also the detection 
of the photo-excited phonon states, describing the probing process by a fully quantum approach, providing in 
this way a direct comparison with the experimentally measured photon number quantum fluctuations of the 
scattered probe pulses [37]. 

In this framework the first step is to adopt a quantized description for the mode-locked pulsed laser fields [38]. 
Each mode of frequency ujj =■ ujq + jS, where cjq is the pulse central frequency, 5 is a constant depending on the 
laser repetition rate and j is an integer, is quantized and described by single mode creation and annihilation 
operators oj and Oj. In this framework ISRS can be modelled by means of an effective impulsive interaction 
Hamiltonian which is descriptive of both the pumping and the probing processes. In both processes two optical 
fields with orthogonal polarizations (denoted with subscript x or y) are involved: two pump fields in the pumping 
process and the probe and the emitted field in the probing process. The interaction Hamiltonian has the form 


'H — + &xj &ljj' b) -f gj ji byj' (&^) + b,xj a}yji & )j 

3,3' = -J 


( 1 ) 


where 2 J-|-1 is the total number of modes within a mode-locked optical pulse, b and fit are the phonon annihilation 
and creation operators, /ij and fis are coupling constants and the functions gj j, take into account the relations 
between the frequencies of the involved fields. 


93,3' 


1 if/=i + f ^^12 

0 elsewhere, ’ 


with fl the phonon frequency. A complete interaction Hamiltonian should contain also second order processes 
involving phonons with opposite momenta. However, since the probe detects only the k ~ 0 optical transition, 
we can make use of an effective Hamiltonian that accounts only for this kind of process. 

The whole theoretical description of the experiment can be rationalized in a four step process as sketched in 
Fig. 3: (i) generation of phonon states in the pumping process, (ii) evolution of the produced vibrational state, 
(hi) probing process and (iv) read out of the emitted photon observables. 

(i) Initially, the sample is in thermal equilibrium and it is described by a thermal phonon state pg, at 
inverse temperature /3. The laser pump pulse is described by a multimode coherent state of high intensity 
\v) = \v-j) ® ® where \vj) are single mode coherent states associated with all the frequency components 

within the pulse. Each \vj) is an eigenstate of the annihilation operator cLj of photons in the mode of frequency 
LOj, ttj \vj) = Vj \vj)- We indicate with v the vector whose components are the amplitudes Vj. The equilibrium 
(pre-pump) photon-phonon state p = |7') {v\ ® pp is instantaneously transformed into p'' =IA pW by means of 
the unitary operator lA = exp{—7H}. Since the pumping operator lA acts on a high intensity photon coherent 
state V, we can use the mean field approximation for the photon degrees of freedom and replace a with v and 
with 13* for both pump modes involved in equation (1), thus replacing lA by 


W = exp|-i [gly pd{v*xj Vyf h' + i3*j, b) + g^y PsiKj ^yj' }• (2) 

3,3' = -J 

The evolution operator generates coherent and squeezed phonon states, respectively, through the linear and 
quadratic terms in the phonon operators b and The initial state p^ contains information about both photons 
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and phonons. Tracing over the photon degrees of freedom, the resulting state pJ'j describes the excited phonons 
brought out of equilibrium by the impulsive pump process. 

(ii) The time evolution of the excited phonons is described by using an open quantum systems approach, 
namely by means of a suitable master equation of Lindblad form [39,40] that takes into account, besides 
the quantum unitary evolution, also the dissipative and noisy effects due to the interaction with a thermal 
environment. 

(iii) The incoming probe pulses are in the multimode coherent state |d). The probing process at time t is 
described by the same impulsive unitary operator lA used for the pump. However, in this case we can apply the 
mean held approximation only to the probe photon operators with x polarization, which correspond to a much 
more intense held than those with y polarization. Moreover, since the probe held is much weaker than the pump 
held, the quadratic terms in the interaction Hamiltonian in equation (1) can now be neglected. The resulting 
unitary operator is 


Wa- = exp{-i||a'||(i(a')&^^ , (3) 

II j=-j 

where a' is the vector with components a' = dy j otxj' and A{a.') is a collective photon annihilation 

operator such that H'l'(a')j = 1 . 

The latter unitary operator acts on a state of the form |d) (d| ® p['j(r). The information about the phonons are 
extracted by measuring the emitted held photons. In particular, the emitted photon state /5i(t) is obtained by 
tracing away the phonon degrees of freedom. 

(iv) The possible quantum features of the phonon state, e.g. squeezing, can be read oh as they are imprinted 
into In particular for each time delay r we can compute the quantities (Ny)^ = (Ht(Q')A(d'))^ and 

-V -V ^2 

A'^Ny = i^y)^ ~ which correspond to the observables measured in the experiment, that are the mean 

value and the variance of the number of photons of the emitted held. The details of the theoretical model are 
reported in (Supplementary Note 3). The theoretical results for Ps = 0 and /is 0 are shown in Fig. 4 together 
with the corresponding wavelet analysis for the variance of the number of emitted photons. The calculations 
reproduce the experimental results in Fig. 2, revealing a 2fl frequency component in the variance, only when the 
pump creates squeezed phonon states (/ig 7 ^ 0). In particular, for /tg 0, the model reproduces the different 
lifetimes between the fl and 2fl components in the variance observed in the experiments. The explicit expressions 
for the theoretically predicted amplitudes of both the frequency components in the variance are reported in 
(Supplementary Note 3), showing that the same damping constant, characterizing the dissipative phonon time 
evolution, contributes differently to the two frequency components giving rise to different decay times. 


Discussion 

The proposed effective interaction model is further validated by a pump fluence dependence study. Fig. 5 shows 
the amplitude of the 2H peak in the Fourier transform of the variance, ATvar, as a function of the pump fluence. 
A fluence dependence study of the fl peak is reported in (Supplementary Fig. 5). The functional behaviour 
obtained from the model predictions (continuous line in Fig. 5) agrees with the experimental data only in 
presence of a pump-induced squeezing of the phonon mode (/ig 0). The increase of the 2fl peak amplitude 
with the pump fluence allows us to give a direct estimation of the uncertainties of the phonon conjugated 
quadratures which are reported in the inset of Fig. 5 for the different excitation fluences (calculation details are 
given in (Supplementary Note 3)). For high pump fluences the uncertainty on one of the phonon quadratures 
falls below the thermal limit at the equilibrium, indicating the squeezed nature of the photo-excited thermal 
vibrational states. Our novel experimental approach allows for the direct measurement of the photon number 
quantum fluctuations of the probing light in the shot-noise regime and our fully quantum model for time domain 
experiments maps the phonon quantum fluctuations into such photon number quantum fluctuations, thereby 
providing an absolute reference for the vibrational quantum noise. The comparison of the predicted noise with 
the experimental photon number quantum uncertainty, measured in shot-noise conditions, allows us to unveil non 
classical vibrational states produced by photo-excitation. A future extension of the model taking into account 
the role of the electronic degrees of freedom would allow to extend such a study from transparent materials to 
complex absorbing systems. 

In conclusion, a Raman active phonon mode has been impulsively excited via ISRS in a a-quartz by means 
of a pump and probe transmittance experiment with single pulse differential acquisition in noise conditions 
limited by intrinsic probe photon number fluctuations. A fully quantum mechanical effective model where both 
phonons generation and detection are studied through the same effective coupling Hamiltonian establishes a 
direct connection between the measured photon number quantum fluctuations of the emitted probe field and the 
fluctuations of the atomic positions in a real material. Our novel approach is used here to reveal distinctive 
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quantum properties of vibrational states in matter, in particular the squeezed nature of photo-excited phonon 
states in a-quartz. Finally, we stress that our innovative approach paves the way for future studies addressing 
the role of unconventional vibrational states in complex systems [3,6], and the thermodynamics of vibrational 
states [41,42] possibly in the quantum regime. 
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Histogram plot 



Delay time (ps) 


Figure 1: Time domain transmittance histogram plot. ATi distribution as a function of pump-probe 
delay for a representative scan 1. (a) For each time delay a color plot of the normalized histogram of = 4000 
acquired pulses is shown, (b) Histogram plot of ATi centered at zero. The data shown are obtained with the 
largest pump fluence used in the experiments (25mj cm“^). 
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Figure 2: Time domain transmittance mean and variance, (a) ATmean (blue curve) and ATvar (red curve) 
as a function of the pump-probe time delay. The zero time is the instant in which pump and probe arrive 
simultaneously on the sample. In the inset a zoom of the variance for the first 3 ps is shown, (b) Wavelet analysis 
(Morlet power spectrum) of the variance oscillating part, (c) Fourier transforms of the oscillating parts of mean 
(blue curve) and variance (red curve). In (a) and (c) the left axis is related to the mean while the right axis is 
related to the variance. 
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Figure 3: Sketch of the four steps effective theoretical model. The steps are indicated with roman 
numbers. The details of the theory for each step are reported in the text. The photon and phonon system are 
denoted with I and II, respectively. 
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Figure 4: Model predictions. Theoretical calculations of the mean value and the variance of the number of 
photons of the emitted held. The left axis is related to the mean while the right axis is related to the variance. 
Comparison between the case with squeezing coupling constant Hs = 0 (a) and ^ 0 (b). A wavelet analysis 
(Morlet power spectrum) of the variance is reported in the inset for both cases. 



Figure 5: Fourier transform 2f2 peak amplitude of the variance. Amplitude of the 2fl peak of the Fourier 
transform of the time dependent variance ATvar- The error bars indicate the standard deviation over 10 scans. 
Comparison between experiments and theory as a function of the pump huence. In the inset the uncertainties 
for the phonon position and momentum operators, calculated from the model, are plotted for the corresponding 
pump huences. 
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Supplementary Figure 1: Scheme of the experimental setup. The sample is depicted at the origin of 
the coordinates. The polarization configuration of pump, probe and emitted held is indicated on top-left. After 
the interaction with the sample the probe undergoes a polarization selection in order to detect the emitted held 
only. The differential acquisition system is also sketched. 
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Shot-noise characterization 



Probe power (mW) 

Supplementary Figure 2; Shot-noise characterization of the detection apparatus. Variance of 4000 
acquired differential pulses as a function of the probe power. The vertical error bars indicate the standard 
deviation over 10 repeated measurements, the horizontal error bars indicate the instrumental error of the 
power-meter. 
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Supplementary Figure 3; Noise characterization as a function of the number of acquired pulses. 

Variance, ATyar, in absence of the pump, as a function of the number of acquired successive pulses. The dashed 
line indicates the shot-noise level. The grey area represents the region in which classical slow noise sources 
contributes to the variance. The error bars indicate the standard deviation over 10 repeated measurements. 
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Supplementary Figure 4: Noise characterization as a function of the unbalance. Variance, ATvar, 
in absence of the pump for a fixed number of pulses {N = 800), as a function of the unbalance between the 
transmitted probe pulse and the reference pulse, expressed as acquired mean voltage ATmean- Grey area: range 
of minimal noise; yellow area: unbalance region in which the noise randomly fiuctuates in a range of minimal 
values; green area: unbalance region in which the experiments have been performed. Notice that we use the 
convention of expressing the voltage acquired for every single pulse acquisition as the sum of the voltages digitized 
for 500 points at 1 GSs“^. For reference, the digitized measurement of a differential pulse with a 1 V (—1 V) 
voltage peak corresponds to a value of 300 V (—300 V). The error bars indicate the standard deviation over 10 
repeated measurements. 
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Pump fluence (mJ cm ^) 

Supplementary Figure 5; Amplitude of the Q peak of the Fourier transform. Fourier transform 
peak of the time dependent mean, ATmean, (blue points) and of the time dependent variance, ATvar, (red points). 
The dashed lines are a guide for the eyes. The error bars indicate the standard deviation over 10 scans. 
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Supplementary Note 1. Polarization dependent ISRS in quartz 


Impulsive stimulated Raman scattering (ISRS) is a non resonant excitation mechanism of lattice vibrations in 
transparent materials by ultrashort laser pulses. It is a four-wave mixing process due to third order polarization 
effects. The Cartesian components of the third order non-linear polarization are given by 

xUl [S 2 {u^ 2 )]k [^ 3 ( 0 . 3 )]/, ( 1 ) 

jkl 


/q\ 

where Xijki susceptibility tensor, fi is the probe field and S 2 and Ss are the pump fields. The susceptibility 

tensor determines the polarization selection of vibrational modes that can be excited via ISRS [1]. 

In particular, quartz Raman active vibrational modes are 4 totally symmetric modes of symmetry Ai and 8 
doubly degenerate modes of symmetry E (transverse and longitudinal) [2]. 

In our experiment, the sample is a 1 mm thick a-quartz, oriented in order to have the principal symmetry 
axis parallel to the probe propagation direction. A scheme of the chosen experimental geometry is shown in 
Supplementary Figure 1. Both pump and probe come from the same laser source, a 250 kHz mode-locked 
amplihed Ti:Sapphire system. The pulse duration is 80 fs, the fractional change in the probe transmission due 
to the pump is of the order of 5% for a pump fiuence of 25 mj cm“^. The excited phonon state is detected via 
the scattering of the probe pulse which arrives on the sample with time delay t with respect to the pump. The 
transmitted light undergoes a polarization selection through a polarizer positioned after the sample. 

The pump direction is almost collinear with the probe one. Assuming that the involved optical fields propagate 
along the 2 direction, we can limit our analysis to the xy plane. In this case the quartz Raman tensors assume 
the form [3] 


A = 


a 0 
0 a 







( 2 ) 


Following the notation in [4], the susceptibility tensor can be expressed as: 




— Ai j Aki -I- Efj eT, -\- E^j Eh , 


(3) 


where each index can assume the values 1,2 associated to the direction x and y respectively. Thus, the 
(3) 

susceptibility tensor Xijki gives rise to a 4 x 4 block-matrix 


ffa^+c^ 0 \ f 0 c^] 

Vo a^-cy [c^ Oj 

fO cy 0 \ 


\ 

/ 


(4) 


The first two indexes refer to the outer matrix elements and describe the polarization components of the emitted 
field (z index) and of the probe field {j index), while the last two indexes (fc and 1 ) indicate the inner matrix 
elements and describe the polarization components of the two pump helds. 

In particular, we are interested in selecting the excitation of an E symmetry Raman mode. For this purpose 
we use a probe linearly polarized along x and we perform a polarization selection after the sample in order to 
detect the emitted field component orthogonal to the probe (along y). This polarization configuration allows the 
selection of the susceptibility matrix elements associated with the involved E phononic mode. Notice that 
such elements vanish when k = I that is when the two pump helds are both polarized along x or along y. Thus, 
in order to activate the process, we need the two pump helds (two frequency components of the same laser pulse) 
to have orthogonal polarizations. This is possible when the pump pulse is linearly polarized along a direction in 
between x and y. In particular, the efficiency of the ISRS is maximal when the pump polarization is at 45 ° with 
respect to the x axis. This is indeed the conhguration we chose and consequently the matrix elements involved in 
our experiment are X 2 iA;/(fc^() ~ getting an emitted held almost collinear with the unscattered probe photons 
and with polarization orthogonal to the probe one. We conhgure a polarizer after the sample in order to transmit 
the emitted held polarization only. The global polarization conhguration is sketched in Supplementary Figure 1. 

Note that a full extinction of the unscattered probe is experimentally not achievable (polarizer extinction 
rate 10®). The residual probe acts as a local oscillator amplifying the emitted held within the total signal [5]. 



Supplementary Note 2. Single pulse differential acquisition system in shot-noise 
limited regime 

The acquisition system is made of a balanced amplified differential photodetector and a fast digitizer {Spectrum 
M3i.2132-exp) with sampling rate 1 GSs“^. The differential photodetector consists of two Hamamatsu S3883 
Silicon PIN photodiodes with 0.94 quantum efficiency connected in reverse bias and followed by a low-noise charge 
amplifier. The photo-currents generated by the two photodiodes in response to a single optical pulse impinging 
on them (transmitted probe pulse on the first photodiode and reference pulse on the second photodiode) are 
physically subtracted and the resulting charge is amplified using CAEN custom designed electronic components. 
In particular, the preamplifier sensitivity is 5.2 mVfC”'^ with a linear response up to about 2 V of pulse peak 
voltage. Note that we use here the convention of expressing the voltage acquired for every single differential pulse 
acquisition, AT^ , as the sum of the voltages digitized for 500 points at 1 GSs"'^. For reference, the digitized 
measurement of a pulse with a 1 V voltage peak corresponds to a value of 300 V. 

Notice that our acquisition apparatus distinguishes itself by avoiding the lock-in amplification and the 
possible artifacts associated with its use [6]: the single pulse differential acquisition system we adopted allows the 
individual measurement of each single transmitted probe pulse, and gives access to the intrinsic photon-number 
quantum fluctuations. 

In order to distinguish the intrinsic noise (shot-noise) from other contributions, we tested the experimental 
set-up in absence of the pump. The variance ATvar of 4000 differential pulses is measured for different powers of 
the probe. The optical noise, shown in Supplementary Figure 2, is linear with a constant offset representing the 
electronic noise. This behavior is characteristic of the shot-noise regime [7]. We chose a probe power within the 
linearity interval (2.5 mW, 0.2 mj cm“^ on the sample; note that this corresponds roughly to 10® photons per 
pulse scattered in cross polarization on the detector). This provides a reference for the shot noise value of about 
1 which is used to benchmark the noise in time domain experiments and avoid additional noise sources. 

The detector shot-noise linearity test demonstrates that we are sensible to quantum fluctuations of the photon 
number. In particular, the shot-to-electronic noise at the maximum probe intensity in the linear regime is 
approximately 10 dB. 

We further characterized our experimental apparatus by measuring the variance, ATvar, in absence of the 
pump, as a function of the number of acquired pulses. The results are shown in Supplementary Figure 3. We 
observe that the variance increases when the number of acquired successive pulses increases, this means that for 
long acquisitions a slow noise contribution makes higher the measured variance (grey area in Supplementary 
Figure 3). Thus, we chose to acquire N = 4000 successive pulses per step in order to guarantee the statistical 
meaningfulness of the data but at the same time avoiding contributions of classical slow noises. Note that the 
time to time noise in the mean number of photons reported in the main text is larger than the value of the 
variance. This is explained by considering that the fluctuations in the mean values are made large by slow noise 
of classical nature (the differences between the measurements scan by scan). On the contrary the variance is 
dominated by intrinsic fluctuations and, as discussed in the main text, it is calculated for every scan separately 
and then averaged. The results reported in Supplementary Figure 3 show that acquiring up to about 10'^ pulses 
guarantees that classical slow noise contributions are excluded. 

Moreover we measured the obtained variance ATvar for a fixed number of pulses {N = 800) in absence of the 
pump as a function of the unbalance between the transmitted probe pulse and the reference pulse, that is as a 
function of the acquired mean voltage ATmean- The results of such characterization measurement are shown in 
Supplementary Figure 4. One can notice that the noise randomly fluctuates in a range of minimal values (grey 
area) for small positive or negative unbalance (yellow area). For larger unbalance the noise starts to increase due 
to artifacts in the amplification process. All the time domain experiments reported here have been performed 
within the region of small detector unbalancing (green area) in order to be sure of working in shot noise limited 
conditions. 

The results of our time resolved experiments show the presence of a 20, frequency component only in the 
variance of the probe photon number. A pump fluence dependent study of the of the amplitude of the 20 peak 
in the Fourier transform of the variance is reported in the main text. For completeness we report here the 
amplitude of the 0 peak in the Fourier Transform of the mean, ATjnean, and of the variance, as a function 

of the pump fluence. The data are show in Supplementary Figure 5. 
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Supplementary Note 3. An effective fully quantum mechanical model for ISRS 


The effective fully quantum mechanical approach to ISRS followed in the main text consists in the pump process, 
the subsequent dissipative, irreversible phonon dynamics and the probe process, all of them being described by 
quantum dynamical maps [8]. 

Before being hit by the pump laser beam described by photons in a multi-mode coherent state \ v) {v\, the state 
of the relevant phonon mode at frequency H is appropriately taken to be a thermal state at inverse temperature 

P 

pp = {l-e-P^) . (5) 

The pump process is characterized by a photon-phonon interaction Hamiltonian of the form 




[djj 




,6t + 


CL'P o CL^ 


^3 “W 


,b) + 


9j,r 9. 


(C- 






( 6 ) 


where /id and /tg are coupling constants, 2J -|- 1 is the total number of modes within a mode-locked optical pulse, 
and the functions take into account the relations between the frequencies of the involved fields. 


I 

93,3' = 


if j' =j + 

elsewhere. 


m 

s 


’ = 1,2 


As stated in the main text, it should be noted that while the linear term involves only the creation of a phonon 
in a single mode at null momentum k, the quartic term are not limited to k = 0 and one should integrate over 
the entire optical phonon dispersion including processes where the momentum conservation is guaranteed by the 
creation of optical phonons with opposite momenta [9,10]. In our effective Hamiltonian we include only a single 
phonon mode at k = 0. This assumption is made in view of the fact that, in the performed experiments, the 
probing process is limited to the linear regime so that phonons at k 0 will not affect the observed photon 
number fluctuations in this configuration. 

Initially, the sample is in thermal equilibrium and it is described by a thermal phonon state pp, at inverse 
temperature /?. The Hamiltonian in (6) generates an impulsive change of the initial photon-phonon state 
\v) {v\ ® pp given by 

p'' =U {\u) {v\ ® pp)U'' = \v) {v\ (g) WppUl , (7) 

where, because of the high intensity of the pump laser beam, we have adopted the mean field approximation and 
substituted the photon annihilation and creation operators by the scalar amplitudes v and i3* and replaced U 
with 


Ujy = exp{-i[ci + c* S + C 2 + C 2 ^^1} (8) 

j 

Cl = /id 9),j' Kj '^V3' (9) 

J 

C2 = /is glj' 13*3 Vyy . (10) 

3,3'= -J 

The pump thus prepares the relevant phonon degree of freedom in a state /jJ) which is obtained from p"' by 
tracing over the photon degrees of freedom: 


p^^, = Tn{p’^)=WppUl 


( 11 ) 


where H and I refer to the phonon and photon system, respectively. 

The linear contribution in the phonon operators is responsible for the displacement of b and while the 
quadratic one accounts for their multiplication by hyperbolic functions and thus for the possible squeezing of the 
corresponding quadratures [llj: 


m 



-b 


1 


(S-l) 



/ cosh(2jc2l) —sinh( 2|c2|)\ 

^_e-d<i>+f) sinh(2|c2|) cosh(2|c2|) / ’ 


( 12 ) 

(13) 


where C2 = |c2|e*'^. In order to write the squeezing matrix S in the standard formalism [11], we can define for 
convenience a complex squeezing parameter 

.7 

i = where r = 2jc2| = 2 |/is| ^ g'jyvlyVyj/, and 3l3 = (j)+-. (14) 

3,3' = -J 
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Notice that the squeezing parameter amplitude r depends linearly on the intensity of the pump pulse and on the 
squeezing coupling constant /ig which weights the non linear term in the interaction Hamiltonian and models the 
material properties involved in the process. 

_ S+fct 


The variance of the quadrature operator B = with respect to the state pJ'j is given by 


A^.B = Tr„ 


- (^Trii (po-B)) = ^ coth [cosh(2r) - sinh(2r) cosi/'] 


(15) 


Then, ior ip = 0 and r large enough, one can make A?c B smaller than 1 /2 which is the shot noise variance of B 
with respect to the vacuum state |0) such that 6|0) = 0. 

The photoexcited phonon state pfj then undergoes a dissipative dynamics that effectively takes into account 
the interaction of the phonons with their environment until, after a delay time t, the target is hit by the probe 
laser beam. The phonon dynamics is considered to be that of an open quantum system in weak interaction with 
a large heat bath that will eventually drive the time-evolving phonon density matrix Pn(t) = Pb(^) to a thermal 
state ppi at temperature T' larger than that of the pre-pump phonon state: P' < p. Such a relaxation process is 
described by a master equation [12,13] for the phonon density matrix ph{t) of the form dtPh{t) = L[/3b(t)j, where 
the generator of the time evolution is given by 


L[/3b(t)] = pb{t)j 

+ A(l-l-n') (^6pb(t)&^-/3b(t)|) 

-I- Xn'(b^ pb{t)b-, pb{t)^^ , 


(16) 


where n' = 


{> n = 


while A is a coupling constant sufficiently small so that the non-negligible 


g/3'O _ I — eP^-\> 

presence of the environment can nonetheless be accounted for, in the so-called weak-coupling limit regime [12], 
by a master equation of the above type. 

The first term of L generates the rotation in time of the phonon mode phase at its own eigenfrequency. The 
second two contributions consist of a so-called noise term h ph{t) bp respectively b^ ph{t) b that has the property of 
transforming pure states into mixed states and of a dissipative term — 116^ S, Pb (t) |, respectively “ 51 bbP , pb (t) |. 
These terms counterbalance the noise by keeping the trace of the time-evolving state Pb(t), and thus the overall 
probability, constant in time. The anti-commutators can be incorporated into the Hamiltonian as anti-Hermitian 
contributions responsible for exponential time relaxation. The structure of L is such that the generated time- 
evolution maps, formally 7^ = exp(tL), compose as a forward-in-time semigroup: 7t ° 7s = 7s ° 7t = 7t+s for all 

> 0. Moreover, Pb(t) = 7t[pb] can be explicitly computed for any initial phonon state pb; all initial states are 
eventually driven to a unique invariant state satisfying L[pb] = 0 that is given by the thermal state pp'. 

Finally, the probe process is again described by the Hamiltonian in equation (6). However, the corresponding 
impulsive unitary operator U = exp {—iH) now acts on a photon-phonon state of the form |a) (d| ® Pn('r). Here, 
d) (d| is the multi-mode coherent state associated with the probe laser beam which contains x and y polarized 
components and is much less intense than the pump one, while Pn(T) is the phonon state dissipatively evolved 
up to the delay time t between pump and probe. Differently from the pump process, the lower probe intensity 
allows one to neglect in T-L the quartic terms responsible for the squeezing effects. Moreover, we can apply the 
mean field approximation only to the field operators with x polarization, since these probe components are much 


more intense than those polarized along y. Then, by replacing Uxj and 
described by 


^ j by a, 

IA-, = exp{—i||d'|| (H(d') fef-I-(d') 6 )} , 
where ^(d') is the collective photon annihilation operator 


xj and apj the probe process is 


(17) 


\\a' 


- E K)' 

j=-J 


a, = Pd 


j'=-j 


Sj',j ' 


(18) 


\ayj) is the 


Then, the probe process affects an initial state |dy) (dyj ® p'ppr), where jdy) = \ay-j) ® 
coherent state involving only the y polarization components such that Uyj jdj,) = Oyj |dy). 

Notice that, unlike in (7), acts on the photon-phonon state as a whole and transforms it into 

Wa' I«S/) (dyl ® Pii{t-)uI, . (19) 

This allows for the quantum features of the phonon state and of its dynamics to be transcribed onto the emitted 
photon state 

Pi{t) =Trii (U^, |dj;)(dy| (» p['i(T)Wt j . (20) 
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Unlike in the semi-classical theoretical approaches to pump and probe experiments attempted so far, one can 
here confront the experimental data not only with the scattered probe beam intensity, namely with the mean 
photon number (Ny)^, where Ny = {a')A(a') and 

(iV,)^=Tr(iV,pi(r)) , (21) 


/s /S -'2 

but also with its variance A'^Ny = i^y)^ ~ {^y)T- 
Then one uses that 

Ul.NyUw = ^+(d')^(a')cos2(||a'||)-f6t6sin2(||a'||)-F Jsin(2||d'||)(A(d')&^+^^(d')&), (22) 

where, given the experimental conditions effectively described by the model, it is plausible to set all amplitudes 
axj = O.X and (Xyj = ay = \ay\ex]3{i6y), in which case 

a' = MdOx = iMdOrrle*® , and ||d'|| = \/if|^da^| , (23) 

where K = 2J + 1 is the total number of modes within a mode-locked optical pulse. 

By denoting with ly = K\ay\'^ the pulse intensity for the y polarization and using that 

i(d') Id,) = Id,), (24) 

one explicitly computes: 

{Ny)x = IyCOs\\\a'\\)+sin\\\a'\\)(b^b)^+^-^ySin(2\\a'\\){e-^^<>'-^yHb^^ (25) 

where {6)^ = Tru 6^ is the expectation value of any phonon operator O with respect to the phonon 

state p\i{t). 

Despite its complicated expression, we report also the number variance A^^Ny predicted by the model, as (Ny)^ 
and A^iV, are the quantities computed numerically in the main text and compared with the experimental data: 

^rNy = I,cos^(||d'||)-f sin^dld'll) (^((6+6)2)^ _ -f sin^dld'H) cos^dld'H) (b^b)^ 

- /,sin2(||d'||)cos2(||dd|) (((St)2)^ _ (^^^2)^ - (6)^)] 

-f /,sin2(||d'||)cos2(||d'||) ( 2 (b^b)^ + l- 2 (b^)^(b)^'^ 

-I- sm(||d'||)cos^(||d'||) (6^)^ - (b)^^ 

+ i^y sin^dld'll) cos(lld'll) [2e-‘(«'-^«) (((6^)^6). - {b^b)^ {b^)^ + ^ (fet)^ j 

- 2e‘(‘^'-^y)(^(b^b%-(brbl(bl + ^(bl^ . 

( 26 ) 


The phononic correlation functions involving b and b^ contribute with oscillations at frequency D while those 
involving r and 6^2 gjyg j-jgg 20 oscillations. Collecting the corresponding coefficients one finds the following 
amplitude for the 20 oscillating components: 


|^2o(r)| = Mi±Mg-Ar Sin2(2||dd|) sinh(2r) , (27) 

o 

where the amplitude of the squeezing parameter r = 2|c2| = 27i'|/is| |j/|2 is obtained from (14) and (10) by putting 
all pump amplitudes equal to v. Moreover we take A to comply with the observed oscillation time-scale and the 
time T > 0 such that At 1. 

In the last figure in the main text we have shown a fit of the experimental results for different pump intensities 
with the functional behaviour of j42n predicted by the model in equation (27). We found an optimal value of the 
coupling parameter /is for which the model agrees with the experiments. We used such a value for computing 
the amplitude r of the squeezing parameter (defined in equation (10)) for all the experimental pump fluences. 
In particular |ii|2 is the number of photons per unit cell per pulse. We then computed the uncertainties in the 
position and momentum phonon operators as in equation (15). The results reported in the main text unveil 
photo-excited thermal squeezed vibrational states. 
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For completeness we also report the explicit time evolution of A'^Ny in terms of both the amplitude of 
the 2 fl frequency component and the amplitude Aq, of the fundamental frequency component; 

AlNy = Ao{t) + An{T)e^^ + (28) 

where the explicit expression for |j4s-i(T)| is; 

kin(T)| = |2e“'^^ sm(2||Q:'|[) sin^(||d'||) [(1 + n' + n — (1 + 2n') cosh(2r)) e*®’’ cosh(r) (29) 

+ (n' — n + (1 + 2 n) cosh(2r)) sinh(r)] 

— sin(2||d'||) (l + 2nsm^(||Q:'||)) (e*®’’cosh(r) — sinh(r)) I , 

where \z\ is the corresponding photo-exited displacement in the phonon and its phase. 

We stress that, if the pump pulse does not generate squeezed phonons, vanishing squeezing parameter (r —> 0), 
the amplitudes of the two frequency components become, 

^o(t)| = |sin(2||d'||) -I- 2sm(2||d'||) sin^(lld'll) (n - (n - n')e“^'^)| , (30) 

|^2n(T)| = 0, 

indicating the absence of the 212 frequency component in the variance in absence of phonon squeezing. 

From equations (29) and (27) one can notice that the damping constant A, characterizing the dissipative 
phonon time evolution between the excitation and the probing process, contributes differently to j42o and to Aq, 
giving rise to different decay times for the two components and reproducing the experimental results. 
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